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Abstract 

The paper deals with a class of cooperative functional differential equations (FDEs) with infinite 
delay, for which sufficient conditions for persistence and permanence are established. Here, the 
persistence refers to all solutions with initial conditions that are positive, continuous and bounded. 
The present method applies to a very broad class of abstract systems of FDEs with infinite delay, 
both antonomous and non-autonomous, which include many important models used in mathematical 
biology. Moreover, the hypotheses imposed are in general very easy to check. The results are 
illustrated with some selected examples. 
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1. Introduction 

This paper is concerned with the persistence and permanence for a class of fnnctional differential 
equations (FDEs) with infinite delay written in abstract form as 

x'i{t) = Fi{t,xt) - Xi{t)Gi{t,Xt), i = l,...,n (1.1) 

where F = {Fi ,..., F„), G = {Gi ,..., Gn) : D IR" are continuous on D C IR x H, and B is an 
adequate Banach space of continuous functions defined on (—oo,0] with values in IR”. As usual, Xt 
denotes the entire past history of the system up to time t, or, in other words, Xt{s) = x{t + s) for 
s < 0. Special emphasis will be given to the autonomous version of (1.1), 

x[{t) = Fi{xt) - Xi{t)Gi{xt), i = l,...,n (1.2) 
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where : rj C -B ^ IR are continuous for all i. 

We shall consider (1.1) as a general framework for some models from mathematical biology, 
therefore only positive (or non-negative) solutions of (1.1) are meaningful. Initial conditions at a 
fixed time a will be taken to be x(s) = ip{s), s < a, with ip : (—oo,0] —>■ IR” continuous, bounded 
and positive (or non-negative). 

The family of FDEs (1.1) is required to satisfy the so-called quasimonotone condition [30] - 
here, we abuse the standard terminology, and these systems will be also called cooperative systems 
(see [30] for a rigorous definition). The main technique employed here is the theory of monotone 
dynamical systems for retarded FDEs well-established in [29,30,34]: the quasimonotone condition 
implies that the solution operator is monotone relative to the initial condition function ip (but not 
necessarily strongly monotone) and allows comparison of solutions between related FDEs. However, 
it is not always clear how to apply the tools of this theory to the situation of infinite delays. In 
particular, for (1.1) it is not obvious which conditions should be imposed in order to guarantee the 
persistence of all solutions with positive, bounded initial conditions. 

From earlier works (see e.g. [16,25,26]), it has been apparent that a rigorous treatment of FDEs 
with infinite delay requires a particular attention to its abstract formulation. In Section 2, we choose 
an admissible phase space B to deal with (1.1), in the sense that it should satisfy some fundamental 
axioms introduced by Hale and Kato [17]. To address the question of permanence, we shall need 
bounded positive orbits to have compact closure in B, so, in addition, we choose a state space that 
is a fading memory space [14,18,27]. The contents and organization of the remainder of paper are 
described below. 

Section 3 is the core of the paper, where the main theorems are established, referring to persis¬ 
tence (Theorem 3.1) and permanence (Theorem 3.2) of autonomous cooperative systems (1.2). A 
key point of our research is that quite reasonable sufficient conditions are enough to guarantee the 
persistence (or even the uniform persistence) of system (1.2) in the entire set of solutions with initial 
conditions that are bounded, continuous and positive. In Section 4, we use the results and tech¬ 
niques in Section 3 to further deduce some persistence and permanence criteria for non-autonomous 
cooperative systems (1.1). We emphasize that, both in Section 3 and 4, no global asymptotical 
constancy is invoked to establish the persistence or permanence of a system. 

Section 5 is dedicated to applications and is divided into three subsections, each dealing with a 
selected important example from population dynamics. In the first example, sufficient conditions are 
established for the persistence and permanence of a cooperative Lotka-Volterra model with infinite 
delay and patch structure, improving significantly some results in a former paper of the author [10]. 
In fact, the study of this Lotka-Volterra model was a main source of motivation for the present study, 
since there was an open problem concerning the persistence left unsolved in [10]. For alternative 
techniques for competitive Lotka-Volterra system, see e.g. [1,8,24,32,33]. The second example 
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refers to an FDE modelling the growth of a single or multiple species divided into n classes and 
following a modihed delayed logistic law; the system has unbounded discrete time-varying delays, 
also includes dispersal terms among the classes, and may be interpreted as a generalization of the 
modified scalar delayed logistic equation proposed by Arino et al. [6]. The last selected example 
concerns an FDE system for the interaction of two species structured into mature and immature 
classes, with the two adult populations in competition; this system is based on Aiello and Freedman’s 
model for a single species [2], and the hnite delay case was studied by Al-Omari and Gourley [5], 
among others. Although the system is competitive, the method in [5] relies on the construction 
of sequences of auxiliary cooperative systems, which are used to prove the global attractivity of 
a positive equilibrium. Here, we show that Al-Omari and Gourley’s method can be extended to 
the case of infinite distributed delay. These examples have been considered by many authors, and 
additional references will be given in Section 5. Many other examples from the literature could be 
analyzed, but the main purpose of Section 5 is to illustrate the application of our main results. 


2. Preliminaries: phase space and notations 


In this preliminary section, we set an abstract framework to deal with FDEs with inhnite 
delay. In view of the unbounded delays, the phase space B shonld satisfy some fundamental axioms 
which gnarantee that the classical resnlts of existence, uniqueness, continuation, and continuous 
dependence of solutions are valid - a subject well establish in the literatnre. Secondly, for our 
purposes, it is desirable that positive orbits of bounded solutions (with bounded derivatives) are 
precompact in B. A convenient choice of B is set below. 

Let g he a function satisfying the following properties: 


(gl) g : (- 00 , 0 ] ^ 

(g2) lim ^ = 1 uniformly on 

“^ 0 - g[s) 

(g3) gis 


[l,oo) is a nonincreasing continuous function, 5 ^( 0 ) = 1; 
oo,0]; 


oo as s ^ —oo. 

For re G N, define the Banach space B = C^, 


= G°(1R") :=U€ C((-oo,0];]Rn 


lim = 0 

s-^-oo gi^s) 


with the norm 



MM 

9{s) ’ 


where | • | is any chosen norm in IR”. To hx terminology we suppose that IR"" is equipped with 
the supremum norm, |x| = jxjoo = max |xi|, for x = (xi,... ,Xn) G IR"^. Consider also the space 

l<i<n 

BC = HC’(1R"') of bounded continuous functions 4> '■ (—oOiO] ^ IR”. Here, BC is considered as a 
subspace of some space (7°, so BC is endowed with the norm of C^. The more explicit notation 
BCg will be also used. 
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The space is an admissible phase space [17,18] for n-dimensional FDEs with infinite delay 
written in the abstract form 

x'it) = ( 2 . 1 ) 

where / : T) C IR x ]R” is continuous and, as usual, segments of solutions in the phase 

space Cg are denoted by Xt- Xt{s) = x{t + s),s < 0. When / is regular enough, the initial value 
problem is well-posed, in the sense that for each (u, (/?) G D there exists a unique solution x{t) of 
the problem x'{t) = f{t,Xt),Xcr = v?, defined on a maximal interval of existence. This solution 
will be denoted by x{t;a,ip) in IR” or xt{cr,ip) in Cg. For autonomous systems x'{t) = f{xt) 
and If G (7°, the solntion of x'{t) = f{xt), xq = if is often simply denoted by x{t;ip) G IR"' and 
Xt{ip) G Cg. When considering more than one FDE of the form x'{t) = f{t,xt) or x'{t) = f{xt), 
the notations x{t] a, ip, f) or x{t] ip, f) will be also used, where the argument / is included to avoid 
any confusion over which EDE is being considered. The space Cg has important properties, namely 
it is a fading memory space - although not always a uniform fading memory space. Eor definitions 
and resnlts, see [14,18,27]. For FDEs x'{t) = f{xt) in fading memory spaces, positive orbits of 
solutions x{t) = x{t‘,ip) with bounded initial conditions (i.e., for ip G BC) which are bounded and 
have bounded derivatives on [ 0 , oo) are precompact in C® - a property that will be used to prove 
our results of permanence. 

Alternatively, one may consider the phase space B = UCg, where UCg = UCg{Si^) :={(/) G 
C’((—oo,0];IR"") : sup^<o < oo, is uniformly continuons on (—oo,0]}, with the norm || • ||g 

defined above. The space UCg is an admissible phase space, in the sense of Hale and Kato [17], 
but it is not necessarily a fading memory space - contrary to what is stated in [19, p. 47]. In fact. 
Haddock and Hornor [14] completely described the functions g satisfying (gl)-(g2) for which UCg is 
a fading memory space; in particular, g must satisfy the condition < g{s) < e~^^^ for s < —M, 

for some 7 i, 72 ,M > 0. For example, (gl)-(g3) hold for g{s) = 1 — s, but g does not satisfy this 
latter condition. If g{s) = e~'^^,s < 0, for some 7 > 0, the spaces Cg and UCg are uniform fading 
memory spaces, and are often denoted by C^ and UCj, respectively. Some of the resnlts established 
here require having precompactness of bounded positive orbits, therefore there is a clear advantage 
in considering the phase space Cg, rather than the nsual UCg. 

We now set some notation. A vector c in IR"" is said to be positive if all its components 
are positive, and we write c > 0. A function ip : (— 00 ,0] ^ IR*^ is said to be positive, with 
notation 99 > 0, if the vectors ip{s) are positive for all s < 0. We define and denote in a similar 
way non-negative vectors and non-negative functions. In the space (7° (or UCj), a vector c is 
identified with the constant function V’('S) = c for s < 0. The non-negative cones of IR” and BC are 
]R!( = {x G ]R” : X > 0} and = H(7_(_(]R’^) = {(/? G BC : (p{s) > 0 for all s < 0}, respectively. 
Motivated by the applications to mathematical biology, we shall take the set 

BCq = {(/9 G BC : ip{s) > 0 for all s < 0} 
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as the set of admissible initial conditions for (2.1). However, in applications, frequently more general 
orders on B are considered, induced by other cones K C B. We use {x)i or simply Xi to designate 
the component i of a vector x € IR"'. If / is a function with values in IR"", it is understood that fi 
means the ith-component of /, for 1 < i < n. 

3. Main results 

In this section, we prove the main results of the paper, about persistence and permanence for 
a class of autonomous FDEs of the form 


Xi{t) = Fi{xt) - Xi{t)Gi{xt), i = l,...,n, (3.1) 

where F = (Fi,..., F„), G = (Gi,...,G„) : BC — IR"" are continuous. It is understood that 
BC = BCg is taken as a subspace of H = G° {or B = UCg). As mentioned above, we are concerned 
with the solutions x{t) = x{t](p) of (3.1) with initial conditions 

Xq = ip with ip G BCq. (3-2) 


We start with a preliminary lemma. 

Lemma 3.1. Consider F,G : BG IR” continuous. Then: 

(i) if F, G are bounded on bounded sets of BG, and x{t) is a non-continuable solution of (3.1) 
defined on [0,a) with a < oo, then limsupj_^„- l®(^)l = oo; 

(a) if Fi{ip) > 0 for all ip G BGj,. with V9i(0) = 0 (1 < i < n), then the solutions x(f) with initial 
conditions ip G BGj^ satisfy x{t) > 0 for t > 0, whenever they are defined. 

Proof. The statement in (i) is a classical result [17,18]; the proof of (ii) follows from [30,34]. ■ 
System (3.1) reads as 

x'it) = f{xt), (3.3) 

where 

/ = (/i, ■■■Jn), fM) = FM) - <i)i{b)Gi{4>) for (f G BG, 1 < i < n. 

In the remainder of this section, F = (Fi, ..., Fn),G = (Gi, ..., G^) are assumed to be con¬ 
tinuous, bounded on bounded sets of BG, and regular enough so that the initial value problems 
(3.1)-(3.2) have unique solutions - which is the case if / is Lipschitz continuous on bounded sets. 
The main assumptions that will be imposed, either in part or in total, are described below. 

(Al) for (j),ip € BG.^-,(j) < if and 4>i{0) = then (i) Fi(4>) < Fi('ip) and (ii) Gi{(j)) > 

Gi(if), i = l,...,n; 

(A2) for each (5 > 0 and each i = 1,.. . ,n, there exists e > 0 such that for all </>, V' G BGj,., (f < if 
with jjV’ — 4>\\g < £ aiid <fi(0) = ifi(0), then Gi{(f) — Gi{if) < 5', 
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(A3) there exists a constant c > 0 such that for ()>,'0 G BC+,(I) < 0 and 000) = 000), then 
F00) - F00) > -c||0 - 010, i = 1 ,..., n; 

(A4) F is sublinear in IR0, i.e., for x G IR0 and a G (0,1), F{ax) > aF{x); 

(A5) there exists a vector v G IR0 such that F{v) — Bv > 0, where B = diag (Gi(0),..., G„(0)); 
(A6) there exists a vector q = (gi,..., g„) G IR0 such that F0g) — qiGi{Lq) < 0 for L > 1, i = 
1 ,... ,n. 

Some comments about the choice of the above hypotheses are in order. 

It is apparent that if the pair (A, G) satisfies one or more of the above hypotheses, the same 
happens for the pair {F,G), where i0(0) = i0(0) + G0O)00O),G00) = Gi{(j)) — G0O), 1 < i < n. 
Hence, without loss of generality, we can take G(0) = 0, in which case B = 0 for the matrix B 
in (A5). Assumption (Al) asserts that F and —G, and thus / in (3.3) as well, satisfy Smith’s 
quasimonotone condition given by (see [30]) 

(Q) for 0,0 G HC+,0 < 0 and 000) = 000), then /00) < /00), i = 1,... ,n. 

This condition implies that the semiflow for (3.3) is monotone on HC+. Here we abuse the termi¬ 
nology, and also refer to systems satisfying (Q) as cooperative systems. Instead of (A2), we could 
simply demand G to be uniformly continuous on BC^, but this requirement is too strong for our 
purposes. Another advantage of the above formulation is that (A2) is trivially satisfied if each G00) 
depends only on 000). It is clear that (A3) and (A4) are satisfied if F is linear bounded; also, if 
F is Lipschitz continuous on BC^, then (A3) holds. Condition (A6) is true whenever either G = 0 
and F{q) < 0, or lim Gi{Lq) = oo, i = 1,..., n, for some vector q > 0. A careful analysis of the 

L—yoo 

proofs below shows that, as alternatives to the sets of constraints (A4),(A5) and (A4),(A6), one can 
simply impose the hypotheses 

(A5’) there is u = (ui,... , Vn) G IR0 such that Fi{lv) — lviGi{lv) > 0 for 0</<l,l<i<n, 
and 

(A6’) there is g = (gi,..., qn) G IR0 such that Fi{Lq) — LqiGi{Lq) < 0 for L > 1,1 < i < n, 

respectively, without the additional requirement of the sublinearity of . Although not essential, 

the sublinearity condition (A4) allied to the quasimonotone condition (Al) implies, however, that 
solutions with non-negative initial conditions remain non-negative. In fact, for (p G with 

Pi{0) = 0, from (Al) we have fi{ip) = Fi{ip) > i0(O), 1 < i < n. On the other hand, (A4) implies 
F(0) > 0. From Lemma 3.1(ii), it follows that x{t; (/?) > 0 for t > 0 in its interval of definition. 

For the definitions of persistence and permanence given below, see e.g. [19,31]. 

Definition 3.1. A system x'{t) = f{xt) with S C BG^ as set of admissible initial conditions 
is said to be persistent (in S) if any solution x{t; ip) with initial condition G S is defined and 


6 



bounded below from zero on [0, oo), i.e., 

liminf Xi{t;ip)>0, 1 < i < n, 

t^OO 

for any if G S. The system is said to be permanent (in S) if it dissipative and uniformly persistent; 
in other words, all solutions G S, are defined on [0,oo), and there are positive constants 

m, M such that, given any ip G S, there exists to = to{ip) for which 

m < Xi{t] If) < M, 1 < i < n, t > to- 

Here, unless stated otherwise, we take S = BCq as the set of admissible initial conditions. 

As observed above, from (A4) we get F{0) > 0. The situation Ti(0) > 0 for one or more 
components of the vector T(0) is not excluded from our setting, although T(0) = 0 is a quite natural 
condition: for many applications of interest, (3.1) corresponds to a population dynamics model, for 
which zero should be a steady state. If F{0) > 0, the persistence of (3.1) follows immediately from the 
quasimonotone condition (Q): with /(O) = F{0) > 0, the solution x^{t) := x(t;0) is nondecreasing 
in t > 0 [30, p. 82] and strictly increasing on an interval [0,e] (e > 0). From the order-preserving 
semiflow, this shows that liminft_,.oo a;i(t; </^) > liminft_>.oo (t) > x^{e) > 0,1 < i < n, for any 
If G BC+. 

In the case of bounded delays, there is an extensive literature using the theory of monotone 
dynamical system to study the persistence of both autonomous and non-autonomous delayed pop¬ 
ulation models. The situation is much more complex for unbounded delays, since some of the usual 
methods do not work unless additional conditions are imposed: in spite of a careful choice of an 
appropriate phase space, as the one set in Section 2, it is a rather difficult task to deal with solu¬ 
tions with initial conditions ip = {pi ,... , pn) in the full set BCq, due to the fact that one may have 
min inf pi{s) = 0. An alternative way to by-pass this difficulty, is to consider only solutions with 

l<i<n s<0 

initial conditions whose components are bounded below and above on (—oo, 0] by positive constants. 
The challenge here is to obtain sufficient conditions for the persistence of (3.1) in BCq - a situation 
not often addressed in the literature, unless there are additional conditions on F, G which allow to 
relate solutions of (1.1) with solutions of an already known permanent system, or the permanence 
appears as a by-product of the global attractivity of a positive equilibrium. 

The main result about persistence is given below. 

Theorem 3.1. Suppose that the initial value problems (3.1)-(3.2) have unique solutions defined on 
[0, oo). If assumptions (Al)-(A5) are satisfied, system (3.1) is persistent (in BCq). 

Proof. For a positive vector v = (vi, ... ,Vn) as in (A5), fix any small > 0 such that 

F{v) -(B + 5I)v > 0. (3.4) 
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Consider c > 0 as in (A3): for 0, V’ ^ BCj^ with (j) <'4’ and 0i(O) = '0^(0) 


^ < Fi{4) - Fi{(t)) < c\\4 - 4\\g, l<i<n. 


By (Al)-(A2), there exists e > 0 such that for 0, 0 G BC-^,(j) < 0, with 000) = 000) and 
110 - 010 < e, 

0 < G00) — G00) < 5/2 for i = 1 ,..., n. 


Step 1. Let <p G BCq be given. Choose Mq > 0 such that 


2||(^||oo 

g{-Mo) 


and 


CM 

9{-Mo) 


< 


S 

2 


min Vi, 

l<i<n 


(3.5) 


where |u| = |u|oo- Note that Fi{v) — {viGi{v/m) + Svi) Fi{v) — (G0O) + S)vi > 0 as m —>■ oo. Fix 
mo G N such that 

Fi{v) — ViGii — v\ — 5ui > 0, i = 1,... ,n,m > mo, (3.6) 


and 


(/9(s) > —u on [—Mo,0]. 
mo 


As before, write (3.1) in the form (3.3). Assumption (Al) implies that (3.1) is cooperative. 
Together with (3.1), consider the auxiliary FDE 


x{t) = f{xt), 


(3.7) 


where = (//,... , /^), //(0) = /00) - 5000) for 0 G BG+,i = 1,..., n. Next, we denote x{t) = 
x{t-, (p, f) and y'^{t) = x{t] ^v, f^) for any m > mo. 

From (3.6) and the sublinearity of F, for m > mo and i = 1,..., n. 



> 


m 


Fiiv) - ViGi[ —V 

m 



> 0. 


Clearly, system (3.7) is cooperative as well, consequently y0(t) is nondecreasing on [0,oo) (1 < i < n) 
(see Theorem 5.1.1 and Corollary 5.2.2 of [30]). In particular, in0>oy0(t) = ^ > 0. 

Step 2. We now prove the following claim: 


x{t)>y'^{t) for t>0,m>mo- (3.8) 

If the claim fails to be true, there are m > mo and to > 0, i G {1,..., n} such that 

Xj{t)>y]^{t) for t e [0,to),j = 1,... ,n, 

Xi{to) = vTito)- 
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Hence 


0 > x'iito) - {yT)'{h) 

= F,{x,,) - + x,{to)[5 - + G,(y-)] 

= - F^) + F^) - F.iyrj + xM [s - G,(xtJ + G^) - G^) + G^{yZ. 

where we take 

mm{ipi{to +s),^Vi), s < -(Mq + to) 


(3.9) 


(/)(s) {(plj ■ ■ ■ ) 4^n)i '/’i ('S) 


and 


'tp{s) = 


^yTi^o + s), -(Mo + to)<s<0 

(p{-Mo), s<-(Mo + to) 


x{to + s), — {Mq + to) < s ^ 0 

Since Xtg(s) > ^{s),'ip{s) > yt7(®) (“00,0] and ^i{0) = ^i(O) = Xi(to) = vTito), by (Al) we have 

Fi{xto) — FAZ) > 0, Gi{'ip) — Gi{yZ < 0. On the other hand, from (3.5) 


m 

^to 


< 


< 


5vi 


^ ~ m i<i<n g{—{MQ + to)) ~ mg{—Mo) '' 2mc’ 


and 


Xto-V\\g= sup 


Z{s + to) -ip{-Mo)\ 
9{s) 


< 


g{-Mo 


< e. 


This implies that 0 < Fi{y'^) — Fi{^) < Svi/{2m) and \Gi{'ip) — Gi{xtg)\ < (5/2. From (3.9), we 
therefore obtain 

But Xi(to) = yi(to) > ^ from Step 1. The above inequality yields 0 > —<5^ + |a:i(to) > 0, which 
is not possible. This proves the claim (3.8). 

From Steps 1 and 2, 


liminf Xi{t) > liminf yZi^) > t = 1 ,..., n, 

t—>-oo t^oo 777,0 


(3.10) 


which shows the persistence of (3.1) in BGq. 


Remark 3.1. Consider an FDE with bounded delays, written in abstract form as (3.1) in the 
usual phase space G = C([—r, 0]; IR"") with the supremum norm. Let (7+ be the non-negative cone of 
G, and take int(C+) = {ip ^ G : (p{s) > 0 for s G [—r, 0]} as the set of admissible initial conditions. 
In this case, we can compare directly the solutions x{t;(p) and x{t; ^n) of (3.1), for m G N small 
enough so that p{s) > on [—r,0]. In this way, the persistence of (3.1) in int{G^) is trivially 
obtained by assuming (Al), (A4) and (A5) (or, alternatively, (Al) and (A5’)). Since the system is 
autonomous, this leads to the persistence of (3.1) in Cq := {p G C+ : (/9(0) > 0}. 
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Corollary 3.1. Consider F : BC^ —>■ IR"" satisfying (Q), (A3), (A4). If there is a vector v > 0 
such that F{v) > 0, the system 

x'{t) = F{xt) (3.11) 

is persistent in BCq. 

In the case of finite delays, again the hypothesis (A3) is not required in the above corollary. It 
should be stressed that, even in the case of finite delays, this corollary provides a better criterion 
for persistence of cooperative FDEs than the one in [35]. 

To address the permanence of (3.1), we start by establishing the boundedness of all solutions. 

Lemma 3.2. Assume (Al)-(A4) and (A6). Then the solutions of the systems (3.1)-(3.2) are defined 
and bounded on [0,oo). 

Proof. Again, write (3.1) in the form (3.3). From the sublinearity of 1^1 , for g = {qi,..., q^) > 

0 as in (A6) and L > 1 we have F{Lq) < LF{q), thus 

fi{Lq) < L[Fi{q) - qiGi{Lq)) <b, i = l,...,n. (3.12) 

From the quasimonotone condition, Xi{t;Lq, f) is nonincreasing in t € [0, oo) and x{t]ip,f) < 
x(t; Lq, f ) if ip < Lq [30]. This proves that all solutions are defined and bounded on [0, oo). ■ 

The main result about permanence of system (3.1) in the full set BCq is stated below and 
requires a fading memory space as a phase space. 

Theorem 3.2. Let F, G be continuous, bounded on bounded sets of BC, and sufficiently regular so 
that the initial value problems (3.1)-(3.2) have unique solutions. If (Al)-(A6) are satisfied, system 
(3.1) is permanent in BCq. Moreover, there are positive equilibria x*,y* of (3.1) such that any 
solution x(^t) of (3.1) with initial condition xq = p a BCq satisfies 

X* < liminfrc(f) < limsupx(t) < y*. (3.13) 

>oo 

Proof. Let q G IR” be as (A6). Since bounded positive orbits are precompact in B = Cg (or 
B = UG.y for some 7 > 0), then x{t;Lq,f) \ y*{L) as t ^ 00 , where y*{L) is an equilibrium of 
(3.1) [30, p. 82]. 

In a similar way, for v = {vi,..., Vn) > 0 as in (A5) and e G (0,1) sufficiently small, we get 

fi{ev) > e{Fi{v) - ViGi{ev)) > 0, i = l,...,n. (3.14) 

Thus x{t',£v,f) is nondecreasing in f G [0,oo). By Lemma 3.2, all solutions with initial conditions 
in BGq are bounded, hence we conclude that x{t; ev, f) x*(e) as t ^ 00, where x*(£) is a positive 
equilibrium of (3.1). 
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Let ip G BCq be given, and take (5 > 0 such that (3.4) holds. From Theorem 3.1 (see (3.8)), 
we also obtain that there exists mo G N with x{t; ip, f) > x{t; ^v, f^) for t > 0 and m > mo. We 
further prove the following claims. 

Claim 1. There is Sq > 0 such that for 0 < e < eo, 

x{t]ev,f)/^x* as t —)■ oo, (3.15) 

where x* is an equilibrium of (3.1). 

To prove this claim, we adapt some arguments in [30, pp. 62]. For similar ideas, see also 
[29,35]. As shown above, if Sq > 0 is such that fi{ev) > 0 for 0 < e < eo and 1 < i < n, then 
x{t;£v,f) x*(e) as t ^ oo, where x*(e) are equilibria, and x*(e) < x* =: x*(£o)- Moreover, 

since {dx/dt){0;£v,f) = /{ev) > 0, then x{t;£v,f) > x{0;£v,f) = £v for alH > 0 and 0 < e < Sq- 
Clearly, there is C > 0 such that x* < l*v, so the suprema 

(5(e) = sup{5 > 0 : x*(e) > (5u} 

are well defined for e G (0,eo]. Note that (5(e) > e > 0. If 5(e) < eo, then x*(e) > 5(e)u implies 
that x*(e) = x(t;x*(e),/) > x(t;5(e)u,/) > 5(e)u, which contradicts the definition of 5(e). Hence, 
5(e) > eo. We now obtain x*(e) = x(t;x*(e),/) > x{t;£ov,f) x* a.s t ^ oo. Consequently, 
x*(e) = X* for 0 < e < eo. 

Claim 2. There is To > 0 such that for L > Lq, 

x{t; Lq, f) \ y* as t —>■ oo, (3.16) 

where y* is an equilibrium of (3.1). 

We argue as in the proof of Claim 1, so some details are omitted. For Lq > 0 such that Lqq > £ov 
and fi{Lq) < 0 for L > Lq, 1 < z < n, we derive x{t; Lq, f) \ y*{L), with Lq > y*{L) > y* where 
y* := y*{LQ) > x*. For L > Lq, define 

M(L) := inf{M > 0 : y*(L) < Mq} > 0. 

If M{L) > Lq, we get y*(L) < x{t-,M{L)q, f) < M{L)q for t > 0, which is not possible from the 
definition of M{L). Hence M{L) < Lq, and we get y*{L) < x{t; LQq, f) —)■ y* as t —)■ oo. This 
proves that y*{L) = y* for L > Lq. 

Claim 3. For x* = x*(£:o) and y* = y*(Lo) as in Claims 1 and 2, the estimates (3.13) hold; in 
particular, (3.1) is permanent in BCq. 
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To prove the estimates (3.13), fix any (/? € BCq and choose L > Lq such that (p < Lq, where 
Lq is as above. From Claim 2 it follows that limsup^_^^ Xi{t] ip, /) < lim^-yco Xi{t; Lv, f) = y*, i = 
1,..., n, hence (3.1) is dissipative. 

Next, we remark that system (3.7), obtained from (3.3) by replacing Gi((^) by Gi((^) + (5 for 
1 < i < n, also satisfies the assumptions (Al)-(A6) for any 5 > 0 such that (3.4) holds. Proceeding 
as in Claim 1, for <5 > 0 small, we get that there exists eo(^) > 0 such that for 0 < £ < £o(^); 

x{t-,ev,f^) Z'x} as t > oo, 

where > 0 is an equilibrium of (3.7). Now, from (3.8), there exists mo{6) G N such that 
for m > mo(S) we have x{t;p,f) > x{t; ^v, f^) =: With m > l/eo(^), we get that 

\\mm.it^oox{t-,ip, f) > We choose e.g. (5^ = 1/fc for fc € N, and note that increases as 
k increases. By the definition of uiq = mo(<5) in the proof of Theorem 3.1 and the definition of 
£o = £o(<^) in Claim 1, we may take m^ := mo((5fc) ^ oo. Given ki,k 2 € N with /c2 > fci, for m G N 
with m > max(£o(<5fci), So(^fc2))) we get . Thus we deduce that the sequence of vectors 

;= increases at k increases, and is bounded from above by y*, thus there exists x* := lim^ X^. 
Moreover, = {X^,... ,X!^) is an equilibrium of (3.7) with 5 = 1/k, i.e., /i(A^) — ^X^ = 0, and 
liminft_,.oo x(t; (/9,/) > X^. By letting k —>• oo, we obtain the estimate liminit^oo x{t; p, f) > x*, 
where x* is an equilibrium of (3.1). ■ 

Corollary 3.2. Under conditions (Al)-(A6), there exists a positive equilibrium of (3.1). Further¬ 
more, if this positive equilibrium is unique, then it is a global attractor of all solutions with initial 
conditions in BCq. 

Remark 3.2. As a consequence of the above proof, the set [p G BC : Eqv < p < Loq} is 
positively invariant and a global attractor for the semiflow of (3.1). However, if the equilibria x*,y* 
in (3.13) are distinct, this set may contain a very complex dynamics, with more equilibria, periodic 
orbits, or heteroclinic orbits (see e.g. [22,23]). 

The next result is important for applications. 

Corollary 3.3. Consider an FDE with infinite delay of the form 

x[{t) = Fi{xt) - Xi{f)Gi{xi{f)), i = l,...,n (3.17) 

where Fi : BC —IR, : [0,oo) ^ IR are Lipschitz continuous on bounded sets of their domains, 
and Gi{0) = 0, 1 < i < n. Suppose that F satisfies (Q), (A3), (A4), and that there are positive 
vectors v,q such that F{v) > 0,F{q) < 0. Then the solutions of (3.17) with initial conditions 
xq = p € BCq are defined on [0, oo) and (3.17) is permanent in BCq. 

Proof. Since Fi,Gi are Lipschitz continuous on bounded sets, initial value problems have 
unique non-continuable solutions. For (3.17), conditions (Al)(ii) and (A2) are trivially satisfied. 
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and F{q) < 0 for some vector g > 0 implies (A6). The latter also implies that solutions are defined 
on [0,oo). The permanence follows from Theorem 3.2. ■ 

4. Persistence and permanence for non-antonomons systems 

We now consider a class of non-autonomous FDEs with infinite delay given by 

Xi{t) = Fi{t,xt) - Xi{t)Gi{t,Xt), i = (4.1) 

for t > to, where F,G : D C M, x Cg ^ IR” are continuous and bounded on bounded sets of D, 
a € ]R, D = (a, oo) x Q, to > a and is an open set containing BG^. Initial conditions are written 
as Xtg = <-p with ip G BCq- In this section, we always assume that solutions to FDEs with initial 
conditions in BCq are unique and defined on [0,oo). 

Theorem 4.1. (i) Assume that there are continuous functions F\G^ : BG —)■ such that 

F\4>)<F{tA), G{t,4>)<G^{<i>) for {t,4>)eD. (4.2) 

If the pair {F^,G'^) satisfies (Al)-(A5), (4.1) is persistent in BGo- 

(a) Assume that there are continuous functions F\F^,G^ ,G^ : BG —)■ IR*^ such that 

F\<i)) <F{t,(j)) <F^{<i)), G\<j)) <G{t,<j)) <G^{())) for it,(j)) e D. (4.3) 

If (F^,G“) satisfies (Al)-(A5) and (F“,G*) satisfies (Al)-(A6), (4.1) is permanent in BCq. 

Proof. In what follows, and without loss of generality, suppose that to = 0. Write (4.1) as 

x'{t) = f{t,xt), t>0, (4.4) 

where fi{t, 0) = E)(t, 0) — (j)i{0)Gi{t, 0), 1 < i < n. Under (4.2), consider the autonomous system 

x'iit) = F({xt) - Xi{t)Gi{xt) =■■ fiixt), i = l,...,n, (4.5) 

and assume that {F\G^) satisfies (Al)-(A5). Since (Q) holds for f \ for non-negative solutions of 
(4.1) we have x{t;0,(p, f) > x{t;(p,f^) {(p G BCq). Theorem 3.1 guarantees that (4.5) is persistent, 
so (4.1) is persistent as well. Similarly, if (4.3) holds we further compare the solutions of (4.1) with 
the solutions of the cooperative system 

x'iit) = Fl'ixt) - Xi{t)G{{xt) =:/("(xt), i = l,...,n, (4.6) 

Since x{t; (p, /*) < x(f; 0, Lp, /) < x(f; ip, /“) for p G BCq, (4.1) is permanent. ■ 

Under conditions (4.3), in fact we obtain explicit lower and upper uniform bounds for the 
solutions x{t;0,p, f),p G BGo, of (4.1): 

X*'’' < lim inf x{t] 0, p, f ) < lim sup x{t; 0, p, /) < y*’“, 

t^-oo 

where are equilibria of (4.5), (4.6), respectively. 

Theorem 4.1 is not easily applicable to FDEs with unbounded time-dependent discrete delays, 
in which case it is better to deal directly with the non-autonomous system. 
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Theorem 4.2. For (4.1), let the pair (F, G) satisfy the assumptions: 

(HI) for (t, (/)), (t, € D with 0 < 4> < then (i) Fi{t,4>) < Fi{t,ip) and (ii) 

Gi(t,(/)) > Giit,-!)), i = 

(H2) for each S > 0 and each i = 1,... ,n, there exists e > 0 such that for all t > to,4>,i(’ G 
BGj^,(f) < i) with ||'0 — (()\\g < £ and (j)i{0) = V’i(O), then Gi{t,4>) — Gi{t,'ip) < 5; 

(H3) there exists a constant c > 0 such that for t > to,(j),'il) ^ BC+,(I) < ip and (pi{0) = V’i(O), 
then Fi{t, p) - Fi{t, ip) > -c\\ip - p\\g, i = l,...,n. 

In addition, suppose that the functions x >-)■ F{t,x) =: F(x) and x ^ G(t,x) =: G{x),x G M”, 
do not depend on t. Then, 

(i) if the pair (F, G) satishes (A4) and (A5), system (4.1) is persistent in BGq. 

(ii) if the pair {F,G) satishes (A4) and (A6), all solutions of (4.1) with initial conditions in 
BGq are bounded. 

Proof. For non-autonomous systems x'{t) = f{t,xt) with f : D <Z Fix G^ ^ IR"" continuous, 
the quasimonotone condition (Q) should be replaced by 

(Q’) for (t, p), {t, p) £ D with 0 < p < p and pi{0) = pi{0), then fi{t,p) < fi{t,p), i = 1,..., re¬ 


condition (Q’) implies that the solution operator xt{cr,p) for x'{t) = f{t,xt) is monotone relative 
to the variable p and that comparison results as in Chapter 5 of [30] are valid. 

(i) As before, write (4.1) as (4.4) and take e.g. to = 0. Assuming (H1)-(H3) means that 
(Al)-(A3) are satisfied with F(-) and G(-) replaced by F{t,-) and G{t,-), respectively (t > 0). In 
particular, (HI) implies (Q’)- 

We now adapt the proof of Theorem 3.1, by replacing equation (3.1) by (4.1) and (3.7) by 

x'(t) = f{t,xt), (4.7) 


where ff{t,p) = fi{t,p) — 5pi{b),l < i < n, and (5 > 0 is such that F{v) — (B + 5I)v > 0 for 
u = (ui,..., Vn) > 0 and B = diag ((5i(0),..., Gn(0)) as in (A5). Clearly, (4.7) also satisfies (Q’)- 
Let ip G BGq be given. It is apparent that the arguments in Step 2 of the proof of Theorem 
3.1 apply with very little change to the non-autonomous case, thus there exists mo G N such that 
x{t) > y'^{t) for t > 0, rre > mo, where x{t) := x{t; 0, ip, /) and y™'(f) := x{t; 0, ^,f^) ■ 

To further prove that inft>o uPit) > 0 for all i and m large, we replace Step 1 of the proof of 
Theorem 3.1 by the following argument: for t > 0 and p > Ay with pi{0) = ^Vi, from (HI) and 
(A4) we have 


ff{t,p) > Fi{t,v/m) - —Vi(Gi(t,v/m) -b > — \Fi{v) - Vi[Gi{v/m) -b <5) 

mV / m L 


> 0 


for m G N sufficiently large and 1 < i < re; from Remark 5.2.1 of Smith’s monograph [30], we derive 
that y™'(t) > Ay for t > 0. This ends the proof of the assertion (i). 
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(ii) In a similar way, under the hypotheses in (ii), a few changes to the proof of Lemma 3.2 
allow us to conclude that all solutions are bounded. ■ 

The permanence result in Theorem 3.2 is not easily adapted to deal directly with non-autono- 
mous systems, since its proof uses cj-limit sets for autonomous FDEs. In fact, an argument often 
used in the proof is that a bounded solution x{t) with monotone components should converge to 
some X* as t —>■ oo, where x* is an equilibrium of the system - which obviously need not happen 
for a non-autonomous system. However, in concrete applications, one might further derive the 
permanence of the system, with explicit upper and lower bounds, as shown in the next section. 


5. Applications to population models 

The above results on persistence and permanence apply to a very broad class of abstract 
FDEs with infinite delay, which include many important models used in population dynamics. An 
advantage of our method is that, in general, the validity of hypotheses (Al)-(A6) and (H1)-(H3) is 
very easy to check. Here, we illustrate the results with some selected examples. 

5.1. A cooperative Lotka-Volterra model with patch structure 

Consider the following cooperative Lotka-Volterra system with patch structure: 


( nOO N 

Pi - fiiXiit) + '^aij / Xj{t - s) dr]ij{s) 
i=i -^0 

/*CXD 

+ '^dij / Xj{t - s) duij{s), i = l,...,n 
Jo 


(5.1) 


where G IR,> 0,dij > 0, r]ij,Uij : [0,oo) —> IR are bounded, nondecreasing functions with 

total variation one, for all i,j. For results on persistence, permanence and stability for autonomous 
or non-autonomous Lotka-Volterra with infinite delays, as well as for a biological explanation of the 
coefficients involved, see [8,10,11,20,21,24,32,33] and references therein. System (5.1) is written in 
the form (3.1) with F), Gi : BC —> IR linear operators given by 


Fpp) = PipipS)+ '^dij / (pj{-s)di^ij{s), 
i=i -^0 

^ poo 

GPP) = iXipip)) - ^ ttij / pj{-s) drjij{s). 
i=i -^0 


(5.2) 


We now insert (5.1) in a phase space C*®, to justify the use of Theorems 3.1 and 3.2. As shown in 
[12,15], given any prescribed d > 0, there is a function g satisfying (gl)-(g3) and such that 


g{-s) dvij{s) < 1 -h (5, 


g{-s)dr]ij{s) <\ +5, i,j = l, 


, n. 
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In this way, Fi, Gi are well-defined by (5.2) as bounded linear operators in (7° [8,18], whose operator 
norms satisfy the estimates ||Fi|| < |/3i| -|- (1 5) dij, ||Gi|| < \fii\ (1 -h (5) hJ = 

1,... ,n. Clearly, assumptions (Al)-(A4) are fulfilled. Now define the matrices 


M = diag (/3i, [dij], N = diag (;Ui,... ,^n) - [dij]. 


For any v G IR*^, F{v) = Mv, G{v) = Nv. For the present setting, we derive the following criteria: 


Theorem 5.1. With the above notations, 

(i) If there exists a positive vector v such that Mv > 0, system (5.1) is persistent in BGq. 

(a) If there exist positive vectors v and q such that Mv > 0 and Nq > 0, system (5.1) is 
permanent in BGq and it has a positive equilibrium. 

(in) In the latter case, if there is a positive equilibrium x* for which Mx* > 0, then x* is a 
global attractor for all solutions of (5.1) with initial conditions in BGq. 


Proof. If Mv > 0 for some positive v G IR*^, then (A5) is satisfied; similarly, Nq > 0 for some 
q > 0 implies (A6). The assertions (i) and (ii) follow from Theorems 3.1 and 3.2. 

We now further assume that Mx* > 0 for some positive equilibrium. In particular, this means 
that Mv > 0 for u = ex* and £ > 0; moreover, for e small, say e G (0,1), as in (3.15) we have 
x{t; ex*) yP X* as t ^ oo. Then any solution x{t) = x{t, (p) {ip G BGq) of (5.1) satisfies the estimates 
(3.13), where y* is also an equilibrium. We need to show that x* = y*. 

In order to simplify the arguments, we consider the ODE system associated with (5.1), obtained 
by taking all the delays equal to zero: 


( n \ n 

Pi - PiXi{f) + '^aijXj{f) j +y^^dijXj{t) 
0=1 / j = l 


i = 1,... ,n. 


(5.3) 


Since the equilibria of (5.1) and (5.3) are the same, we only need to show that solutions of (5.3) satisfy 
limsupj_,.j^ x(t) < X*. For a positive solution x{t) of (5.3), define Lj = limsup4_,.oQ(xj(t)/xp, and 
choose i such that Li = maxj Lj. Choose a sequence —>• oo with xptk) —^ LiX* and x'^tk) —^ 0 
as k ^ oo. Using (5.3), this yields 


0 < LiX* [ Pi - piLiX* -f Lj 

i=i 


O'ijX j 


+ Li d, 
i=i 


tjXj 


= Li{PiX*i p'^dijX*) + Llx*{-piX* p'^aijX*) 
i=i i=i 

n 

= Lpi- Li){P,x*i + y] d,jX*) = Lpi - Li){Mx*)i. 
i=i 


Under the condition Mx* > 0, this leads to Li < 1. From (3.13) and the definition of Li, we get 
X* < liminft^oo 3:j(t) < limsupj_,.j^ xj (t) < x* for 1 < j < n. ■ 
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Remark 5.1. As expected, Theorem 5.1 shows that the global dynamics of (5.1) depend 
heavily on the algebraic properties of the matrices M,N. Condition Nq > 0 for some vector q > 0 
is equivalent to saying that is a non-singular M-matrix; the dissipativeness follows from this 
property. When M is an irreducible matrix, there is a positive vector v such that Mu > 0 if and 
only if s{M) > 0, where s{M) is the spectral bound of M, dehned by s{M) = max{Re A : A G a{M)}. 
If all the coefficients jdi are positive, clearly Mx* > 0 if x* > 0 - so, if > 0 for some q > 0, 
a positive equilibrium of (5.1) must be a global attractor. Theorem 5.1 signihcantly improves the 
results in [10] in what concerns the persistence of cooperative Lotka-Volterra models (5.1): here, 
the persistence is established in BCq, a question which was left as an open problem in [10]. On 
the other hand, [10] deals with questions of extinction, persistence, global asymptotic stability for 
Lotka-Volterra systems (5.1) without any prescribed signs for the coefficients I3i,aij. If there is at 
least one negative coefficient Oij, then (5.1) is not cooperative, and the theory of monotone systems 
cannot be applied directly. In this scenario, auxiliary cooperative FDEs can be used, and results 
of comparison with cooperative systems invoked, to draw conclusions about the non-cooperative 
system - a technique often exploited in [10]. This technique will be illustrated in Subsection 5.3. 

Consider now the following non-autonomous version of the Lotka-Volterra system (5.1): 


/ ^ poo 

x[{t) = Xi{t)(l3i{t) - iii{t)xi{t)+ 'y'aij{t) / Xj{t-s)dr]ij{s) 

V Jo 


1*00 

+ 'y^dij{t) / Xj{t — s) di'ij{s), t > 0,i = 1,... ,n 

i=i -^0 


(5.4) 


where f3i{t),aij{t),dij{t) are continuous and bounded on [0,oo), aij{t),dij{t) are non-negative 
and r]ij, Uij : [0, oo) —)■ IR are bounded, nondecreasing functions with total variation one, for all i,j. 
A straightforward application of Theorems 4.1 and 5.1 gives the following criterion: 


Theorem 5.2. Under the above conditions, define the n x n constant matrices 


= diag (^^,..., PJ + d^j , M“ = diag {(3^,. + 


Hj 


=diag ■ ■ ■, hj - 




where we use the notations 

/ = inf/(^), 7 = sup/(t), (5.5) 

~ *>0 t >0 

for a function f : [0, oo) —>■ IR bounded. Then: (i) if there exists a vector v > 0 such that M^v > 0, 

system (5.4) is persistent in BCq; (ii) if there exist vectors v > 0 and q > 0 such that M“u > 0 and 

N^q > 0, system (5.4) is permanent in BCq. 
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5.2. A system of modified logistic equations with unbounded delays 


Consider the classical delayed logistic equation, also known as Wright’s equation, given by 
x'{t) = rx{t){l — x{t — t)/K), where x{t) is the population of a species at time t, r the intrinsic 
growth rate of the species, K the carrying capacity, and r the maturation delay. In Arino et al. [6], 
the model 


x'{t) 


7/Ux(t — r) 

— l)x(t — r) 


fix{t) — Kx‘^{t), 


(5.6) 


where 'y,n,K,T > 0, was derived as an alternative, and more realistic, formulation for the delayed 
logistic law. The coefficients in the logistic equation are related to the ones in (5.6) by r = 7 — 
(for 7, the birth and mortality rates, respectively) and K = (7 — |u)/k. More recently, Bastinec et 
al. [7] proposed a more general non-autonomous model with multiple time dependent delays: 


x\t) = 


ak{t)x{t - Tfc(t)) 
^ 1 + /3k{t)x{t - Tk{t)) 


— Il{t)x{t) ~ K{t)x^{t), t>0, 


(5.7) 


where ak,n '■ [0, 00) ^ (0, 00), /3fc, ;U, Tk ■ [0, 00) ^ [0, r] (for some r > 0) are continuous, 1 < k < m. 
The question of the permanence of (5.7) was raised (but not studied) in [7], and a criterion for it 
given by the author in [9], with explicit uniform lower and upper bounds for all positive solutions. 

Here, we pursue the research, and generalize the scalar equation (5.7): not only unbounded 
time-varying delays are allowed, but also we consider n classes (of one or multiple species) following 
the ‘modified delayed logistic equation’ (5.7), with dispersal terms among the classes. This leads to 
the system 


(^ik{t)Xi{t — Tik{t)) ^ J (4. _ l4\\ 

1 < R., (f\\ +2^dij{t)xj{t aij{t)) 


1 + /3ikit)Xi{t - Tik{t)) 

- Hi{t)xi{t) - Hi{t)x‘l{t), t > 0, i = 1 ,..., n 


(5.8) 


where all the coefficients aik (t), [dik(t) , dij(t) , , Hi(t) and delays Tik{t), aij(t) are non-negative and 

continuous functions in t € [0, 00), fc = 1 ,..., mi,i,j = 1 ,..., n; the time-dependent delay functions 
Tik{t),o'ij{t) are possibly unbounded. The migration rates of populations moving from class j to 
class i are given by dij{t), with aij{t) the time-delays during dispersion; the instantaneons loss 
term —dji{t)xj{t) may be incorporated in the term —fij{t)xj{t) of the jth-equation. For biological 
reasons, we usually consider dii{t) = 0; in any case for each i G {1 ,..., n} the term dii{t)xi{t — aii{t)) 
may be indued in the first snm on the right-hand side of (5.8). 

We now established sufficient conditions for the persistence of (5.8) in BCq. This example also 
illustrates how to combine the techniques in Theorems 4.1 and 4.2. 

Theorem 5.3. Suppose that aik, ldik,dij, pi, Hi,Tik,<Tij '■ [0,oo) ^ [0,oo) are continuous functions, 
•= ORkit) and Ki{t) are bounded below and above by positive constants, and I3ik,dij,pi 


18 



are bounded, k = 1,..., = 1,... ,n. If there exists a positive vector v = {vi,... , Vn) such 

that 

Hv > 0, (5.9) 

where H is the n x n matrix H = diag {A^ — 7I1, • • • ,^4^ — -|- \_dij\ for 

(*) ’ Ai = inf Ai (t), Pi = sup pi {t), 

t>0 t>0 t>0 

then all solution of (5.8) with initial conditions in BCq are bounded below and above by positive 
constants. 

Proof. As in (5.5), we set the notations / = inft>o/(t), / = supj>o/(f) for a function / : 
[0,oo) ^ [0,oo) bounded. Condition (5.9) reads as 

n 

{Ai-lIi)vi + '^di^Vj > b, i = l,...,n. (5.10) 

i=i 

Effecting a scaling of the variables, Xi{t) = Xi{t)/vi, we obtain a new system with the form (5.8), 
where I5ik{t),dij{t), Ki{t) are replaced by = Vil3ik{t),dij{t) = Vjdij{t)/vi, ki{t) = ViKi{t) for all 

i,k,j, and all other coefficiens are the same. Dropping the hats for the sake of simplification, we 
can therefore assume that (5.8) satisfies condition (5.10) with v = (!,...,!). 

Define rik{t,x) := > 0,x > 0, and note that 0 < ■^rik{t,x) < olik- Write (5.8) as 

x'i{f) = Fi{t,xt) - Xi{t)Gi{t,Xi{t)) =: fi{t,xt), t > 0,i = 1,... ,n, 


with 

rrii n 

Fi{t,(j)) = + '^dij{t)(j)j{-aij{t)), t > 0,0 G BC, 

k=l j=l 

Gi{t,x) = Pi{t) + Ki{t)x, t > 0,x G IR. 

Define also the functions F\t,())), F^{t,(l)) and G*(x),G“(x) whose components are given by 


Fr{t,cf) = j2j 


= y it)), 

'^ik (^)) 


+ y dijfji-aijit)), t > 0, 0 G BG, 


Clearly, 


G\{x) = p. + KiX, Gf{x)=Jii+'KiX, X G IR. 


Fl{t,cf) < F,{t,(f) < i0“(0 0), G\{x) < G,{t,x) < Gf{x) 


for t > 0,0 G BG+,x G 1R+. On the other hand, choose g satisfying (gl)-(g3), and insert (5.8) is a 
space Gg. The functions F\t,(l)),F'^{t,(l)) are uniformly (for t G [0,oo)) Lipschitz continuous and 
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nondecreasing with respect to (p in the non-negative cone and G[{x), G“(x) are increasing for 

X > 0. Therefore the pairs (F*, G“) and (F“, G*) satisfy (H1)-(H3). Note also that x ^ F\t,x),x ^ 
F'^{t,x) are autonomous for x G IR*^, and that F*(x) := x),F“(x) := F'^{t,x) are sublinear. 

Next, consider the auxiliary systems 


x'iit) = F-{t,xt) - Xi{t)Gi{xi{t)) =: fi{t,xt), i = (5.11) 

x'iit) = F^{t,xt) - Xi{t)G\{xi{t)) =: f^{t,xt), i = (5.12) 

By results of comparison of solutions, 

x{t; 0, ^p, f) < x(t; 0, ip, f) < x(t; 0, p, /“), (5.13) 


for p G BGq. From (5.10) with u = (1,..., 1) =: 1 G IR!( and e > 0 small, we have 


rrii n 

F‘{el) - ono) = e -Tt, + Y.i„] >0. i = 


^ 1 + 


n. 


i=i 


Theorem 4.2(i) yields the persistence of (5.11), and thus the persistence of (5.8) as well. Moreover, 
lim2;_^oo G((x) = oo, 1 < i < n, which implies that (F“,G*) satisfies (A6). Theorem 4.2(ii) and 
(5.13) allow us to conclude the boundedness of solutions x{t; 0, p, /) on [0, oo), for any p G BGq. ■ 

With an additional condition on the (possibly unbounded) delay functions, the permanence of 
(5.8) is established, with explicit lower and upper bounds for the asymptotic behavior of solutions. 

Theorem 5.4. Assume the hypotheses of Theorem 5.3, and in addition that lim {t — Tik{t)) = 

t—^CO 

lim {t — crij{t)) = oo for i, j = 1,... ,n, k = 1,... ,mi. Then, system (5.8) is permanent. Moreover, 
for V = [vi,... ,Vn) as in (5.9), the solutions x{t) = x{t;0,p) of (5.8) with p G BGq satisfy the 
uniform estimates 


mo < liminf(xi(t)/ui) < limsup(xi(t)/ui) < Mq, z = 1,..., n, 

t^oo 


(5.14) 


where 


and 


Mq = max lim sup 
l<i<n t->oo 


l<i<n vfKi(t) 


aik{t) - Pi{t) ] +y^^dij{t)i 
/ j=i 


\k=l 


(5.15) 


mo = min lim inf 


1 


l<i<n t^oo V^Ki{t) 


E 


(0 


^ 1 - 1 - l3ik{t)viMo 


- hi(t) ] +'^dij{t)t 


i=i 


(5.16) 


Proof. For x{t) := x(t]0,p, f), set Xj := liminf 4 _,.oo Xj(t), Xj := limsupj_^j^ Xj(t), 1 < j < n. 
From Theorem 5.3, 0 < x^- < Xj < oo for all j. Consider i such that Xi/vi = maxi<j<„(xj/uj), for 
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u > 0 as in (5.9). By the fluctuation lemma, take a sequence {tm) with tm —>■ oo, x[{tm) —>■ 0 and 
Xi{tm) —>■ Xi- For any e > 0 small and m sufficiently large, we have Xi{tm — Tik{tm)) <Xi + £ and 
ViXj{tm — (Jijitm))/vj <Xi + e for all k,j. Recalling that the functions rik{t, ■) are nondecreasing 
on X G [0, oo), for sufficiently large m we derive 


Oi(tm) <(Xi + £)'^j 


^ik 


^ m ){xi+e) 


^ ^ (^m) (^m); 


j = l 




Xi-\- e 
Vihvi(trn) 


Y^Viaik{t 

m ) + ^ ^ dijitr, 
J = 1 


k'i (in 


,k=l 


Ki (R 


~Xi(tm) X^(tm) 


Taking limits m —>■ oo, e —>■ 0^, this estimate yields 


0 < limsup 

t—¥CO 


ViKi{t) 


i=i 


k=l 



thus 

aik(t) - + 

This leads to xjjvj < Mq, for 1 < j < n and Mq as in (5.15). 

To prove the uniform lower bound given by (5.14), (5.16), we reason along the lines above, 
and some details are omitted. Let i G {1,... ,n} be such that Xi/vi = mini<j<„ and take a 

sequence (sm) with Sm —>■ oo, x[(sm) 0 and Xi(sm) —^ x^- For e > 0 small and m large, we get 


Y^dij(t)vj 

i=i 


Xi < lim sup — 

t^oo 


E 

\k = l 


^ x^-e 



^ik 

1 H“ Piki^m^i^i 


1=1 






k'i(^m) 

(^m) 


Xi(Sfn) Xi(Sjji} 


Since Xj < ViMo, taking limits m ^ oo, e —>■ 0+, this leads to 


0 > X j lim inf 


_J_ ^ 


0!ik{Sm)Vi , , \ ^ , / \ 1 


and therefore x^/vi > mo for mo as in (5.16). 


Note that when /dik(t) = 0 in (5.8), the lower bound mo in (5.16) does not depend on Mq. 
In the case n = 1, we get a criterion which improves and generalizes Theorem 3.2 in [9]. 


Corollary 5.1. Consider the scalar equation (5.7), where ak, h, [a, k : [0,oo) ^ [0,oo) are con¬ 
tinuous and bounded, Tk : [0,oo) ^ [0,oo) are continuous, 1 < k < m, and K(t) > Ko,A(t) := 
Ym^=i ^k(i) > Ao for some constants ko,Ao > 0. If 
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then (5.7) is persistent (in BCq) and all positive solutions are bounded. If in addition lim {t — 

t—^OO 

’Tki't)) = oo for k = 1,... ,m, (5.7) is permanent and positive solutions x(t) = x(t; 0, (p) (p G BCq) 
satisfy 


mo < lim inf re (t) < lim sup re (f) < Mq, 

^ t —^oo 


where Mq = 


1 / m \ 1 / 

limsup V ak{t) - p{t) , mo = liminf V - 

t^oo K{tj J K{t) 1 


akit) 


+ (5k{t)Mo 


-hit) 


(5.17) 


5.3. Permanence for a competitive compartmental system with infinite delay 

In 1990, Aiello and Freedman [2] proposed and studied a model for a single species with imma¬ 
ture and mature stages and a discrete time-delay: 

rt'(t) = aumit) - yniit) - ae~^'"ujnit - r) 
u'^it) = - r) - /3u^(t), 

where a,l3,y,T > 0. In (5.17), Ui,Um stand for the immature and mature populations, respectively, 
r is the maturation time since birth, a is the birth rate for the species, and /3,7 the death rates 
for matures and immatures, respectively. Since the second equation is decoupled from the first, to 
describe the qualitative behavior of solutions to (5.17) it is sufficient to study the second equation 
of the system. This model has received great attention over the last decades, a large number of 
generalizations has been derived, and many aspects of their dynamics analyzed, in several contexts. 
Besides the early works [2,3,13], see [4,5,28] and references therein for more results on related models. 
A natural extension is to introduce non-constant delay; in fact, it may even be infinite, leading to 
a second equation given by 

nco 

u'mit) = a / fis)e~^^Umit -s)ds- jdul^it) 

do 

for some / summable with f{s) ds = 1. More recently, models with two or more stage-structured 
species have been proposed. In [5], two species, structured in matures and immatures and in 
competition, were considered. Disregarding the immature populations, this leads to a Lotka-Volterra 
type model of the form 


pCO 

u'i{t) = a\ / — s) ds — I3iu\it) — ci_ui{t)u 2 {t) 

Jo 

pOO 

U^(t) = 02 / f2is)e~'^'^''u2it - s)ds - P2ulit) - C2Ui{t)u2it), 

Jo 


(5.18) 


where aj, l3j,yj,Cj > 0 and the kernels fj : [0,oo) —>■ [0,oo) are continuous (or at least measurable) 
with fjis) ds = 1, j = 1, 2. However, due to technical difficulties, in [4,5] the authors restricted 
their analyses to the case of kernels /i(s), f 2 is) with compact support - of course, if [0, r] contains 
both the supports of fi and /2, the framework is restricted to a system of differential equations with 
finite distributed delay on [—r, 0]. The theorem below was proven by Al-Omari and Gourley [5] for 
the case of a finite delay. 
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Theorem 5.5. [5] Consider (5.18) with aj, /3j,'yj,Cj > 0, fj : [ 0 , 00 ) —>■ [0,oo) continuous, fj{s) = 0 
for s > T for some r G (0, 00 ), and fj{s) ds = 1 , j = 1, 2. Assume that 


C 2 ai [ /i(s)e ^^®ds</ 3 ia 2 / f 2 {s)e ds, 
Jo Jo 

C 1 C (2 [ / 2 (s)e“'^^® ds </ 32 ai [ /i(s)e“'’'^® ds. 
do Jo 


(5.19) 

(5.20) 


Then u* = (u ^, ), where 


(Jih — 


(J1IJ2 — 


- - (l 3 ia 2 [ f 2 {s)e ^^''ds-C 2 ai f /i(s)e ds 

— C 1 C 2 \ do do 

- - (l 32 ai [ fi{s)e~^^''ds - Cia 2 [ f 2 {s)e~^‘^''ds 

— C 1 C 2 \ do do 


(5.21) 


is the unique positive equilibrium of (5.18), and u* is a global attractor of all solutions to (5.18) 
with initial conditions in Cq := {(p : [—r, 0] —> IR^ | p is continuous and (^(0) > 0}. 

Here, as an application of the techniques in Section 3, we prove that the above result is still 
valid for the case of r = 00 . We use the ideas and arguments in [5], inserting them in the present 
framework, which enables us to deal with the infinite delay. Clearly (5.18) is not cooperative, and 
consequently Theorems 3.1 and 3.2 cannot be applied directly. 

Consider (5.18) with kernels /i ,/2 : [0,oo) —>■ [0,oo) continuous with L^-norm equal to 1, and 
assume (5.19)-(5.20) with r = 00 . Choose d > 0 such that 


C2^d/3i + ai j fi{s)e ds^ </3ia2 j f 2 {s)e ds, 

Cl ^d/32 + 02 J f 2 {s)e~^^''ds^ <( 820:1 J fi{s)e~^^''ds. 


(5.22) 


(5.23) 


Fix any solution u{t) = u{t; p) of (5.18) with p G BCq. In (7° or UC-^ with 0 < 7 < min( 7 i, 72 ), 
we write (5.18) in the form (3.1) with n = 2, where F = (^ 1 ,^ 2 ), G = (Gi, G 2 ) are linear functions 
given by 

pOO 

Fi{(j)) = ai /i(s)e“^*'*0i(-s) ds for (/> =((/>!, ^ 2 ) G HG, i = 1, 2, 

Jo 

and G(0) = G((/)i(0), (/) 2 ( 0 )), with 

Gi(u) =/3iui + C 1 U 2 , G 2 (u) =/ 32 W 2 + C 2 U 1 for u = (ui,U 2 ) G IR^. 

Since F(0) > 0 for G ddG_|_, from Lemma 3.1 the solutions u{t) of (5.18) with initial conditions in 
BGq are non-negative in their maximal interval of existence. For such solutions, we derive 


pCO 

u[{t) < ai / fi{s)e~'^^^ui{t — s) ds — (liUi(t) 

Jo 

pCO 

U 2 {t) < 0(2 f 2 {s)e~'^^''u 2 {t — s)ds - (l 2 ul{t), 

Jo 
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and we now compare with solutions to the auxiliary cooperative system 


pCO 

u\{t) = ai / — s) (Is — 

Jo 

/ CO 

/ 2 (s)e“^"% 2 (i - s)ds- /32ul{t). 


(5.24i,J 


System (5.24i^it) satisfies (Al)-(A6), so in particular it is dissipative in BCq. Thus, u{t) satisfies 

limsupui(t) < i = 1,2, 

t—>-oo 


where = (uJ’“,U 2 ’“) is the unique positive equilibrium of (5.24i^„), with coordinates given by 

pOO 

= Ms)e-'^-^ds, i = l,2. 


For t > 0 sufficiently large, we get 


pco 

u[{t) > ai / fi{s)e~^^^ui{t — s) ds —/3iu1{t) — CiUi{t){ul’^ + 5) 

Jo 

pOO 

U 2 {t) > a 2 / f 2 {s)e~^'^‘'u 2 {t - s) ds - (d 2 ul{t) - C 2 U 2 {t){u\’'^ + S), 

Jo 


and compare u{t) with solutions to the new auxiliary system 


/ OO 

fi{s)e~^^^ui{t — s) ds — / 3 iUi{t) — + ( 5 ) 

,.oo (5.24i,0 

U2{t) = a2 / f2{s)e~'^‘^''u2{t - s)ds - P2ulit) - C2U2it){ul’'^ + S). 

Jo 

System (5.24ij) has the form (3.1) and satisfies (Al)-(A4) and (A6). Also, for (5.24i^j) the matrix 


B in (A5) reads as B = 


F{u) — Bu 


_ ( Cl(u 2 ’“ + ( 5 ) 


C2(w}’“ + <5) 


. Formulae (5.22), (5.23) yield 


= {oiiui J /i(s)e ds - ci{u 2 '^ + 5)ui,a2U2 j f 2 {s)e ds) - C 2 (uJ’“ + <5)^2^ > 0, 

for any u = (ui,U 2 ) > 0. In particular, (A5) holds. From Theorem 3.2, we now obtain 


liminf Ui{t) > u\'\ f = 1, 2, 

t—^CO 


where = (uj’^u^’*) is the unique positive equilibrium of (5.24i^i), given by 

1,1 ai/o°“/i(s)e“^i®(is-ci(u2’“ + ( 5 ) ^ ^ 02/o°°/2(s)e“^2^ ds) - C2(%’“ + < 5 ) 

“■ =-A-' ““ =-S- 
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Next, we observe that, for t > 0 sufficiently large, the solution u{t) satisfies 


pOO 

u'i{t) < ai / fi{s)e~^^^ui{t — s) ds — — ciui{t){ul’^ — S) 

Jo 

pOO 

U 2 {t)<a 2 / f 2 {s)e~'^'^''u 2 {t - s) ds — j52U2{t) - C 2 U 2 {t){uY - 6), 

Jo 

and compare it with solutions of the second-upper auxiliary system 

pOO 

u'i{t) = ai / ui{t — s) ds — I5iu\{t) — ciUi{t){vJ 2 ’^ — 6) 

Jo 

poo ’ 

U 2 {t) = 02 / f 2 {s)e~^'^''U 2 {t - s) ds - P 2 ul{t) - C 2 U 2 {t){u\’'' - 6) 


(5.24 


2,u) 


Jo 

which satisfies (Al)-(A6). Proceeding in this way, in [5] the arguments above were iterated and 
auxiliary cooperative systems (5.24^^„), (5.24„^i) with positive equilibria respectively, con¬ 

structed, providing upper and lower bounds for u{t): 


< liminf Ui(f) < limsup ni(t) < z = l,2, n G N. 

t^OO t^oo 

Moreover, explicit recursive formulae can be obtained for the sequences showing that u^’’' 

increases and decreases, i = 1,2, with lim^ = lim„u”’“ = u*. See [5] for details. Thus, 

Theorem 5.6. For (5.18), Theorem 5.5 is valid with r = -|-oo and Cq replaced by BCq. 
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